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WALL BOUNDARY EQUATIONS WITH SLIP AND CATALYSIS 


FOR MULTICOMPONENT, NONEQUILIBRIUM GAS FLOWS 

By Carl D. Scott 
Lyndon B. Johnson Space Center 


SUMMARY 


Boundary equations are obtained for a low Reynolds number, high 
enthalpy gas flow where velocity slip and temperature jump are impor- 
tant. The formulation treats a multicomponent gas mixture that moy be 
in nonequilibrium with finite-rate catalytic reactions occurring on the 
wall. A first-order distribution function is used to include multicom- 
ponent diffusion. A simplified gas/wall interaction is assumed where 
individual atoms or molecules either reflect specularly off the wall or 
stick and are fully accommodated. The boundary equations are obtained 
for use in flow-field calculations such as finite -difference time- 
dependent methods or other methods that are appl .cable to low-density 
flow regimes. 


INTRODUCTION 


The interaction between a high enthalpy, low Reynolds number gas 
flow and a relatively cold surface can be significantly influenced by 
the characteristics of the surface. This interaction establishes bound- 
ary properties for the gas flow and can influence the overall flow field. 
Surface characteristics are important, for example, in the high altitude, 
low Reynolds number flight regime for prediction of Space Shuttle orbiter 
entry heading and for interpreting data obtained in ground tests of space- 
craft entry thermal protection systems. 

In low Reynolds number hypersonic flows such as at high altitudes 
or in low-density arc-heated wind tunnels, the boundary condition assump- 
tions that are adequate for the higher Reynolds number regimes may no 
longer be adequate for describing the flow parameters at a wall (solid 
boundary). Specifically, the density at the wall is sufficiently low 
that the velocity and temperature of the wall are no longer the same as 
those of the gas immediately adjacent to the wall. These phenomena are 
referred tc as velocity slip and temperature jump. 
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Figure 1 is a plot of stagnation pressure as a function of total 
enthalpy showing typical entry trajectories for the Space Shuttle and 
Apollo vehicles. Regions where arc-heated wind tunnels typically oper- 
ate for simulating entry heating and the approximate region in which 
slip effects become important are also shown. It can be seen that there 
is a region of the shuttle trajectory in which slip may have some sig- 
nificant effect on the boundary conditions. Moreover, there is a large 
region of a typical arc-tunnel operating map for which slip is significant. 

For a real, chemically reacting gas, it is necessary to formulate 
a boundary description in which the pressuie and species concentration 
also "jump” at the wall because of rarefaction effects and chemical re- 
actions on the surface. This report is concerned specifically with cata- 
lytic recombination reactions on the surface. 

As the density decreases, the continuum flow equations that ade- 
quately describe higher density flow sire no longer adequate close to the 
wall because the mean free path becomes long compared to characteristic 
lengths associated with significant changes in macroscopic flow proper- 
ties. The flow in a region next to the wall having a thickness on the 
order of a mean free path (the Knudsen layer) cannot be described by the 
Navier-Stokes description (ref. l). Gradients near the wall cause the 
velocity distribution function to deviate significantly from equilibrium. 
Rigorously, one should calculate the flow using the Boltzmann equation 
because there are not sufficient collisions between the gas molecules 
for the Navier-Stokes approximation to be valid. The solution to the 
Boltzmann equation in the Knudsen layer with suitable kinetic boundary 
conditions at the wall must then be matched to the solution of the 
Navier-Stokes equations in the bulk flow. The laws governing the inter- 
action of the gas with the surface influence the boundary conditions and 
must be taken into account. This procedure can be done simply through 
the use of accommodation coefficients. 

Without actually solving the Boltzmann equation in the Knudsen layer, 
it is possible to find an approximate Navier-Stokes solution to the flow 
in the Knudsen layer using suitable slip conditions (refs. 2 and 3). 
Patterson (:. ef. 2) obtained a Navier-Stokes type solution to the Boltzmann 
equation that was of the form of the equilibrium velocity distribution 
function times a power series in the three-dimensional velocity space. 


f(v) = f°(v)fl + 


a. v. 
i i 


+ ± 

2 


a v v 

U i J 




a. .» v.v.v, 
ijk i j kJ 


( 1 ) 


where f(v) is the velocity distribution function; f°(v) 


librium velocity distribution function; 


V 


a. 


ij* 


and 


a i,1k 


is the equi- 
are 
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and 


are 


coefficients of the distribution function; and v. , \ , 

1 J K 

velocity vector components. (Summation convention is used here.) 
Patterson obtained the coefficients in terms of gradients of the macro- 
scopic variables. To obtain the slip relations, he assumed a fully ac- 
commodating surface (diffuse reflections) and equated the net fluxes of 
mass, momentum, and energy to the difference between the incident and 
reflected fluxes, respectively. Shidlovskiy (ref. 3) performed essen- 
tially the same type of calculation except that he allowed for a non- 
fully accommodating surface (specular and diffuse reflection). In both 
of these treatments, a single-species gas was assumed. 


Flow-field and boundary-layer calculations have been made that in- 
corporate limited combinations of the boundary conditions considered in 
this report. Treatments of slip effects can be found in references 1 
to 3- Davis (ref. U) included first-order slip effects in temperature 
and velocity for his reacting binary mixture flow past a blunt body. 

This included calculations for catalycity limits of zero and infinity. 
Inger (ref. 5) included a variable catalytic wall-boundary condition 
for stagnation point, nonequilibrium flow of a binary gas past a sphere. 

He assumed the following: (l) a constant Prandtl number, Schmidt number, 

density viscosity product, and specific heat across the shock layer; 

(2) negligible thermal diffusion; (3) a thin shock layer model; and 
(k) no slip at the wall, but slip at the shock wave. Blottner (ref. 6) 
used a finite-difference method to investigate the vjscous shock layer 
at the stagnation point of a blunt body. An air gas model was used with 
finite reaction rates and transport properties based on Yun and Mason 
(ref. 7) and Yos (ref. 8). Blottner did not consider slip, and he assumed 
a fully catalytic body. Blottner reviewed the literature relative to 
blunt-body flow problems (ref. 6). 


No one has included the effects of slip, multicomponent diffusion, 
and wall catalycity in a single-flow- field problem. In this report, 
boundary relations that incorporate all of these effects for application 
to flow-field calculations are derived. Boundary conditions for the 
temperature, velocity, pressure, and species concentrations are obtained 
for a multicomponent mixture at a wall of arbitrary catalycity. To be 
general enough for low-density continuum flows, slip is also considered 
herein. 


SYMBOLS 


A. 

l 



a. 

i 


>a ij ,a ijk 


coefficients of the distribution function for itu 
species 


coefficients of the distribution function 



1 * 


a iO ,a il’ b iO’ c iO 


(J) 


b 

C. 

1 


distribution function coefficients defined in 
reference 9 for ith species 

impact parameter 

mass fraction of species i 



E 

E 

F 

F 

f. 

l 

f? 

i 

g 

I, I 

k 

k 

v 

t 

M 

M 

m 

i 

n 


diffusion coefficient 
thermal diffusion coefficient 

diffusion vector of jth species 

net energy flux from translational energy 
energy flux from translational energy 
net general flux 
general flux 

distribution function of ith species 
equilibrium distribution function of ith species 

initial relative speed during a collision 

component terms in integration over distribution 
function 

Boltzmann constant 

wall catalytic rate constant 

order of surface recombination reaction 
net mass flux 
mass flux 

mass of ith species 
number density 
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P 


net momentum flux 


P 


momentum flux 


P 


pressure 



collision cross section 




factors in the Sonine polynomial expansion per- 
taining to mass fluxes or species defined in 
appendix A 


Q°° Q 01 Q 10 Q 11 ^ 
lj ,y ij ,H 1J ’ im 


factors in the Sonine polynomial expansion per- 
taining to temperature gradient defined in 
appendix A 


Q 00 R 
\}’ iO 


m 


T 

t 


U 

V 


i 


factors in the Sonine polynomial expansion per- 
taining to velocity gradient defined in 
appendix A 

separation distance of atoms during collision 
distance of closest approach 

temperature 

time 

unit tensor 
thermal velocity 


diffusion velocity 

v velocity 

Vq bulk velocity 


W 


dimensionless velocity = 



V 
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w 

a 

Y 


Y 



n 

0 


0 1 


y 

p 

$ 


<t>(r) 

♦ (V) 

x(e>'b) 


Affixes : 

*.i.J 

s 

v/ 


rate of mass production in volume eleme> of ith 
species 

energy accommodation coefficient 

dimensionless relative velocity in collision, used 
in appendix A 

recombination coefficient 

Kronecker delta 

viscosity 

fraction of incident atoms that are diffusely 
reflected 

fraction of incident atoms that stick 

reduced mass 

density 

perturbation part of distribution function 
interaction potential 

a general property 
collision deflection angle 

collision integrals 


species indexes 

property at edge of Knudsen layer 
wall property 
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Subscripts : 
a 

i.J.k 

x,y»z 

y 

1 1 ,X,Z 

Superscripts : 
q,£,m,n 

s,£ 

+ 

t 


atomic species 
vector component indexes 
component directions 
normal component 
tangential components 


exponents used in flux integrals defined in 
equation (Cl) 

indexes appearing in collision integral 

incident flux 

specularly reflected flux 

ANALYSIS 


The approach to the slip uoundary condition problem follows closely 
that of Shidlovskiy (ref. 3), although additional factors are included 
herein to account for wall reactions. The slip conditions are taken 
across the Knudsen layer, which is on the order of one mean free path in 
thickness. In. figure 2, the Knudsen layer and the jump in temperature 
T are shown. The jump is treated as discontinuous across an arbitrarily 
thin Knudsen layer, and no account is taken of variations of the velocity 
distribution through the Knudsen layer as would be required by the more 
rigorous approach of Kogan (ref. 1). It is assumed here that the dis- 
tribution function near the wall can be described to first-order accu- 
racy by the so-called Navier-Stokes approximation as used by Shidlovskiy. 
The use in this report of a Chapman-Enskog type distribution function 
for a multicomponent mixture obtained by the variational method of 
Hirschfelder, Curtiss, and Bird (ref. 9), however, is a deviation from 
the procedures of Shidlovskiy. The advantage of this usage lies in the 
fact that the Chapman-Enskog distribution function accounts for diffusion. 

The interaction model does not distinguish between energy accommo- 
dation a and momentum accommodation e. For many materials, a and 6 
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are nearly equal arid, have values close to unity; therefore, only a number 
accommodation is considered. By this, it is assumed that the gas atoms 
(molecules) interact with the wall as follows: each atom either reflects 

specularly (is not accommodated at all) or sticks to the wall (is fully 
accommodated). The fraction 1-1 .at sticks is caxled 0'. Of those atoms 
that stick, a certain number . .act with other gas atoms on the wall (usu- 
ally by recombination). The fraction of incident atoms that reacts (recom- 
bines) is called y. Those that do not react leave the wall fully 
accommodated to the wall temperature. This process will be called diffuse 
reflection m this report. The fraction of incident atoms that is dif- 
fusely reflected is 0. This interaction model is considered to be ade- 
quate because the data required for a more detailed model are very 
limited. The change in internal energy of the reflecting molecules, 
such as vibrational and electronic excitation, is neglected. 


The Distribution Functions 

As stated previously, the velocity distribution functions used here 
are those for multicomponent mixtures 



(V)[ 1 + 


1 



(2) 


where 


*.(V) 



In T - B. :Vv. 

l 0 



(3) 


The Maxwellian velocity distribution for the ith species is f?(V). The 

1 ) 1 

coefficients A^, B^, and u are functions of the dimensionless 

velocity = ym^/(2kT)V^ and may he defined as 


\ - [to * “nd • "i 2 )]"i 

B, = b iW.W . - ~ W/Uj 
i i0\ i i 3i/ 


( 5 ) 
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C. 

1 


(* 


) 



( 6 ) 


where a.-., a.-, b.*, and c.- J are constants determined from the 
iCr il iO* lO 

Hirschf elder, Curtiss, and Bird (ref. 9) variation problem in the first 

approximation for a mixture and U is the unit tensor. These constants 
are functions of the collision integrals and are related to the transport 
properties. See appendix A for the explicit fc^m of these constants. 

The vector d 4 is related to the diffusion velocity of the ,]th 
species in the mixture and is defined as 




V In p 


(T) 


where n. and m are, respectively, the number density and mass of the 
J J 

jth species, n is the total number density n = y~n , , p is the total 

j J 

mpss density p = Jja m , and p is the total pressure p » p . . 

j J J j 0 


The 


•fe w *2h. 

total mass averaged velocity = 1/p nijVj where v 


is the total velocity v j = v o + s P ec i es averaged over 

the distribution function and V is the therma.. velocity. 

J 


The Balance Equations 

To obtain the slip relations at a wall, it is necessary to consider 
the balance of normal fluxes of mass, momentum, cna energy at the wall 
for each species. The net flux equals the difference between the inci- 
dent flax and the outgoing flux at the wall. In the fol?.owing equations, 
the net flux will be denoted by an uppercase block letter and the inci- 
dent and outgoing fluxes by uppercase script letters. An upward arrow 
denotes flux caused by specularly reflected particles , and a downward 
arrow denotes incident fluxes. A superscript w denotes outward fluxes 
of wall-accommodated particles. 
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Consider a property <f>(V) such as mass, momentum, and energy, and 
the normal flux of that property across the surface 


F = 



4>(V)f s fv) d 3 


V 


( 8 ) 


is the net flux of A Cartesian coordinate system is used, and 

the uirection normal to the surface is y. The distribution function 

fit) is that just outside the wall on the outer surface of the Knudsen 
s 

layer. 


00 Q CO 


F * = j j Jv y 4>(V)f*(v) d^ 


(9) 


— CD —CO —00 


is the incident flux. 


00 00 03 


= J j J V y 4>(v)f*fv) d^V 


( 10 ) 


—oo Q —oo 


is proportional to the specularly reflected flux. 


00 00 oo 



<f>(V)f w (V) d\ 


( 11 ) 


is proportional to the diffusely reflected flux. These fluxes are shown 

schematically in figure 2. Hie distribution function f vV) is that 

of a gas at equilibrium with the wall temperature. Hie distribution 

function f*(V ,V ,V ) = f*(V ,-V ,V ) because the sign of V 
s x y* z s x y i y 
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changes upon reflection. If $(V) is the mass, then the normal mass- 
flux balance at the vail is: net mass flux equals incident mass flux, 
leaving mass flux. For atomic species 


M i * H i* * t 1 - e 0 M i + ( e ; - W 


(12) 


Because it is assumed atoms are being consumed at the wall by catalytic 
recombination, the net mass flux ^ 0. The first term on the right 

is the mass flux of incident particles, the second term is the flux of 
particles specularly reflected, and the last term is the mass flux of 
particles diffusely reflected (those that stick minus tho^e that 
recombine ) . 

Relations similar to equation (12) for the fluxes normal to the 
surface of the normal component of momentum, the tangential component 
of momentum, and the energy, respectively, are obtained as follows. 


IX * IX * IX 

i i i 

IX i * IX i * 5X 

i i i 

DrS'Mf- 

i i i 






+ 


ix -xr 



i 



(33) 

(14) 

(15) 


By evaluating the integrals represented by equations (9), (10), 
and (11), it can be shown that M* = -M*, P*| | = -P*| | , E* = -E*, 

P I . w = 0, and P* = P f . Equations (12) to (15) then can be reduced to 
ill iy iy 


M. 

1 


9 ; M i * 



(16) 



00|CVJ 
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Du - e o p ir + Z ( 9 i ■ ^ V = pS (11) 


Z p iii * Z e i p ii 


(18) 


Eh - 5 X * 1 > - ^ * 


(19) 


The flux balance equations can v expressed in terms of the flux 
integrals to obtain the following, in which k is the Boltzmann constant. 


Mass of species i: 


00 00 00 


—(Irff f f ♦•.<*>] A 

' ’ *,00 »0Q «cn 


2 oo Q oo 


= 0 




, m n W /2kT \ 2 r r r 

/f/'V 

/-\ * m n m 


-w 2 3 

e W d w 


•OO 0 — 00 


( 20 ) 
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Normal momentum: 




W y 2 e _w2 ^l + *_(W)J d\ 


W 2 e _w2 d 3 W (21) 

y 


Tangential momentums: 





( 22 ) 
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Energy: 



(23) 


w 

The wall number density n. is defined following the convention of 

Shidlovskiy (ref. 3) such that the factor (0 1 - y) appears in the last 
terms of the preceding equations. 


The Slip Equations 

When the integrals in equations (20) to (23) are evaluated and 
after the terms have been rearranged, the following relations for the 
slip quantities are obtained. All accommodation coefficients 6! are 
equal to 0 ' . 


ro|u> 
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Concentration slip (single species): 



(2k) 


Pressure slip: 


^ * - i * *«)», * .-p’ - £*».' 

rz^ 


'iO 


Jy 


( 25 ) 


The viscosity n = x kT^n-b... 

c ^ X lU 

Velocity slip: 



(26) 
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The net energy flux E is defined in appendix B. See appendix C for 
the integrals necessary for determining fluxes. 


The Relation Between Catalytic Recombination 
Coefficient y and Rate Constant 

Using the results obtained for concentration slip, an unresolved 
conflict in the literature may be resolved. The relation between the 
recombination coefficient y and the rate constant k^ will be ob- 
tained. Two expressions have frequently appeared in the literature for 
first-order recombination. One is 


k 


w 



( 29 ) 



IT 


(e.g., refs. 10 and ll), and the other is 

fkT 

“v * < 3 °> 


(e.g., ref. 12). The purported difference (ref. 10) between these two 
expressions stems from inclusion of the diffusion velocity in the deri- 
vation of equation (30). 

These relations can be obtained by using the flux and slip equa- 
tions given previously in this report. The net mass flux of atoms to 
the surface is given by 


M 


i 



(31) 


This is equated to the expression defining the catalytic rate constant 
k^. The rate of consumption of atoms at the wall from surface recombi- 
nation is 


M. 



(32) 


where £ is the order of the surface recombination reaction. For a 
first-order reaction £ - 1, therefore 


M. 


, w 
= -k n, m. 
w i i 


(33) 


The minus sign indicates that the flux is in the direction opposite to 
the outward surface normal. Substituting the expression for the mass 
fluxes (eqs. (Bl) and (B 2)) into equation (31) obtains the relation 


1 3T 

" T 3y a i0 




(3U) 
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Substituting the expression for the net mass flux (eq. (Bl)) into equa- 
tion (32), then 


1 



Substituting equations (34) and (35) and solving for k , then 

w 


1 



If slip and shear are neglected, then this expression reduces to equa- 
tion (30). However, if the species density slip equation (eq. (24)) is 
substituted in equation (36), the expression 


1 



is obtained with no simplifying assumptions. It can be seen that by 

neglecting slip but including diffusion, one obtains the usual form for 

the relation between k and y (eq. (30)) when diffusion is included. 

w 

However, when both diffusion and slip are considered, then the simple 
form (eq. (29)) is obtained. Therefore, the latter form (eq. (29)) is 
to be preferred because it is more general. These expressions imply 
that k never approaches infinity; that is, there is an upper bound 

„n the rate constant which is controlled by the thermal velocity at the 
wall. However, it is sometimes convenient to assume k w ■+ <*> to sim- 
plify the boundary equations for fully catalytic cases. 
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The Species Boundary Equations 

To obtain the boundary condition for species concentrations at the 
wall, one may begin with the species continuity equation. 


TT + , { oc i('o + *i)] 3 “i (38) 

where c^ is the mass fraction of the ith species, t is time, and 

w^ is the ith species volume production term. Although the problems to 

which the boundary conditions will be applied may be time dependent, it 
is assumed that the solutions to be reached are steady state. Then, the 
time- dependent terms of equation (38) go to zero. In addition, on the 

wall, the species volume production term w. vanishes. Therefore, the 
term in brackets is a constant vector. 1 

The normal (y) component is then the net species flux to the surface, 
which has been defined previously. 


pc. s V = p c w k 
i y w i w 


(39) 


The normal bulk velocity v^ is zero at the wall. The preceding is the 

same as the expression obtained in equation (35) for the net species flux 
to the wall. 


S “ w w 

-M. = n. m.V = p ck = n, m.k 
i i x y w i w i i w 


(4o) 


By expressing c 


(J) 


iO 


in terms of the diffusion coefficient D 


ij 


and in terms of the thermal diffusion coefficient D. (ref. 9 - p. 479) > 
then 


3 In 


3y 




i— i- D 

0 ij 


3y \n p / 


3 In p 


ay 


= k n. W m, 
w i i 


( 41 ) 
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If all gradients except the concentration gradient are neglected and if 
a binary gas is assumed, then equation ( Ul ) reduces to 


2 

n m 2 m 1 


2 1^ vn / _ w 

D_ ^ — = k n 

p 12 3y w 1 1 


In terms of the mass fraction of species 1 and neglecting slip, then 




was used here. Equation (U3) is the form of the atom boundary condition 
usually employed in flow-field calculations for binary mixtures. Equa- 
tion (Ul) is the more general boundary condition to be used with slip, 
gradients in T and p, and multicomponent mixtures. 

g 

It may be more convenient to express equation (Ul) in terms of n^ 

and y rather than n. W and k . Equation ( 31 ) is solved to obtain 

1 v 
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DISCUSSION 


From the results obtained in the previous section, the wall boundary- 
conditions can be calculated for a multicompouent gas reacting catalyti- 
cally on the surface, with slip and temperature jump. Equations {2k) 
to (28) and (1+5) form a coupled set of boundary conditions to be used 
with flow-field calculations. However, ? 'y not be necessary to use 
the pressure slip equation. A constant rature boundary condition 

is a natural stipulation for these equat.. nowever, for an adiabatic 

wall boundary equation, the equation for transfer at the wall must 

be set equal to zero, thus adding an additional simultaneous equation ~o 
the set. The wall temperature is then not given as a boundary condition. 
Because the accommodation coefficient 9' is based on the rather crude 
model that either a molecule is not accommodated or is fully accommodated, 
and because for many gas/surface interactions the coefficients are nearly 
equal to unity, it is convenient to assign the accommodation coefficient 
;he value unity. If the gas/solid interaction law is known with suffi- 
cient accuracy, it would be more appropriate to reformulate the balance 
equations taking these known accommodation coefficients into account. 

It may be necessary to define individual momentum accommodation coeffi- 
cients for both tangential and normal momentum accommodation. Rigorously, 
the accommodation coefficients may be functions of velocity and should 
be included in the integrals over the distribution functions . 


To assess the importance of various terms in the general boundary 
equations will require a systematic study in which the flow properties 
and boundary conditions are varied. In many cases, the equations may 
be simplified by choosing appropriate transport property approximations 
or by neglecting various terms in the equations. Coupling of these 
equations with a time-dependent technique such as that of C. P. Li 
(ref. 13) would yield a good calculational procedure to investigate the 
influence of the boundary conditions on the heat transfer and flow. 

Such a procedure should also provide a set of more realistic boundary 
conditions for low Reynolds number hypersonic flows. 
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CONCLUSIONS AND RECOMMENDATIONS 


The equations describing the boundary conditions for a multicompo- 
nent, nonequilibrium gas reacting catalytioally on a surface with tem- 
perature jump and slip are derived on the basis of a simplified 
gas /surface interaction model. These equations form a simultaneous 
set which must be coupled with a flow-field calculation. The equations 
may be used with various calculational techniques in an iterative fashion. 
They are convenient for a time-dependent numerical flow-field calcula- 
tion technique where the time dependence of the boundary properties is 
not of interest. As the flow- field solution approaches steady state, 
the boundary properties approach those given by the boundary constraints 
determined by the flow and the boundary conditions 

To assess the importance of various terms in those general boundary 
equations, a systematic study should be conducted in which the flow 
properties and boundary conditions are varied. Coupling of these equa- 
tions with a time- dependent technique would yield a gooa calculational 
procedure to investigate the influence of the boundary conditions on the 
heat transfer and flow; it should also provide a set of more realistic 
boundary conditions for low Reynolds number hypersonic flows. 


Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, December 18, 1973 
986-15-31-07-72 
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Figure 1.- Plot of stagnat^n pressure as a function of enthalpy 

showing Space Shutf j. 1 c Trajectories with typical arc- 

jex operate - me?, *i;ere slip flow begins for a 

0.3-meter-ra a 0.05-meter-radius sphere are 

denoted. Ike is typical of arc-jet test models. 














F* = incident flux 

F* = specularly reflected flux 

F W = diffusely reflected flux 



► T(y) 


Figure 2.- Diagram of the Knudsen layer showing general fluxes and 
coordinate axes. The temperature as a function of normal distance 
is schematically cverlayed. 
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APPENDIX A 

CONSTANTS APPEARING IN THE DISTRIBUTION FUNCTION 


The general velocity distribution function for one species cf a 
multicomponent gas is given in Hirschf elder , Curtiss, and Bird (ref. 9) 
as 


f t (v) = 


fj(v) 


t* 


^(V) 


] 


(Al) 


In the section entitled "Analysis," it was shown th-.t 



-[ s i0 + *u(f - w i 2 )]“i- 7 ln T 

- b io("i"i - | “i 2 “)= % 




W. »d 


(A2) 


where = ynn/(2kT) V^. This distribution function is % first-order 

approximation sometimes called the Navier- Stokes approximation. The 

coefficients a. ft , a.,,, t. _, and c. ft are found by a variational 
iO* il* 10* lO * 

technique in which they are solutions to sets of simultaneous equations 

(ref. 9).. This appendix expresses the coefficients in terms of solutions 

to these sets of equations. These solutions are expressed in terms of 

Is l) 

the collision integrals 0. v * . 

X J 


a 


ij 


* l) = VUi" / e_Y V S+3 Q (£) (g) dy 


(A3) 


where p is the reduced mass of the colliding particles x and j, 

2 2 

y is the dimensionless relative velocity in collision (y = ug /kT), 
and g is the initial relative speed. 
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oo 



is the collision cross section at collision deflection angle 


x(g,b) = TT - 2b 




dr 


(A5) 


where i|>(r) is the interaction potential, b is the impact parameter, 

and r is the distance of closest approach. The collision integrals 
m 

have been evaluated for a number of interacting ,pecies as a function of 
temperature and are available from several sources (e.g., refs. 8 and 1*0. 

The solutions to the variational problem are from reference 9 where 
simultaneous equations are to be solved. The coefficients of the diffu- 
sion terms in the distribution function axe 



Vjo 


(h) 



(a6) 


where 





* 




n m 

JJ 





(1,1) 
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(h) _ 




(AT) 

(A8) 
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The set of simultaneous equations for the coefficients of the temperature 
gradient terns in the distribution function is 


^TAj^'v = E im B £ {0 ' 1 ' 

J m'=0 


(A9) 


where 


= 8^-=^ — rLm./V. - 6..\ - n.m./l - 6 ..Via., (l,l) 

1J t “O r lj J v J J v /J 


(A10) 




(All) 


qIO = il** 

ij m i \) 


(A12) 


Q 11 

^iJ 




6 ij ” 6 iki 


’(-i ♦ \) 

- ’VV 1 ’ 21 * vV 1 ' 31 


') K 6 ”/ * 5 "k 2 )“lk (1 ' 


1) 


( 4 1J * *Jk) ^j^ik 


( 2 , 2 ) 

(A13) 


and 


*Im * * S ml I 1 “t Jzf' 


(AlU) 
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The set of simultaneous equations for the coefficients of the 
locity gradient terns in the distribution function is 

XXj b J0 = ' R i0 

j 


where 



and 


R 


iO 


-5n. 


ve- 

(A15) 


(A16) 


(ait ) 
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APPENDIX B 

FLUXES IN TERMS OF EVALUATED INTEGRALS OVER 
THE DISTRIBUTION FUNCTION 

The fluxes of a property <f>(V ) and the form of the distribution 

function f(V) are defined and described in the section entitled 
"Analysis." This appendix is a compilation of the forms of these fluxes 
that have been integrated over the distribution functions. 


MASS FLUXES OF SPECIES 


Net 


m.n. 
H i 



9 In T 
3y 


a i0 + 


n ]C c i 


(j) d 

io jy 


(Bl) 


Incident 


, m.n. S /2kT\ 2 b. n /av n 

* _ ill s 1 , . lO I Ox . 

i - - 2 WyT/ 1 y ~ + 

(j) 

d inr n o ^ i n 


3v„_ 9v 


02 ? 

3z 9y 


3y y 


3y “iO 2 fc f #w i0 d jy 
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Specularly Reflected 


(B2) 


m; = - m 


( B3 ) 
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Diffusely Reflected 



NORMAL MOMENTUM FLUX 


(B4) 


Net 


■ 4 ■ ¥(■ 


3v_ 3v_ 

Ox Oz 

3x 3z 


2 !Isl1 

3y 




(B5) 


Incident and Specularly Reflected 
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3z 3y 
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P iy xy p i |^2 6 \ 
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r v Ox 

6 ' 

l 3x 
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Diffusely Reflected 
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TANGENTIAL MOMENTUM FLUX 
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Diffusely Reflected 
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32 


ENERGY FLUX 
Net 



s 


Specularly Reflected 

E* = - e [ (B15) 


Diffusely Reflected 



(Bl6) 
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APPENDIX C 

INTEGRALS NECESSARY FOR DETERMINING FLUXES 


The normal fluxes of mast, momentum, and energy are defined in 
terms of integrals over the velocity distribution function. These inte- 
grals consist of terms that are various velocity moments of the distri- 
bution function for the nonuniform, multicomponent gas and that are 
expressed in terms of the general expressions 


00 CO 00 


■/// 

_co .00 .00 


dW dW dW w q w “w n e~^ 
x y z x y z 


(Cl) 


for the net fluxes, where q, -£, m, and n are exponents hau-ng 
values 0, 1, or 2 and W*" = + W^ 2 + W^ 2 , 


00 0 00 


['*' = f f f dW dV dW w q w “w n e' 

J J J x y z x y z 


n -W 2 


(C2) 


.00 .00 .o 


for the incident fluxes, and 


'-/[I 


2 

dW dW dW wSi “w n e _Vr 
x y z x y z 


(C3) 



3h 


for the specularly reflected fluxes. Specifically, these integrals are 
the following. 


INTEGRALS OF THE ENTIRE VELOCITY SPACE 


(NET QUANTITIES) 
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INTEGRALS OVER THE LOWER HALF SPACE (INCIDENT QUANTITIES) 


III V V ^ - 5 
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INTEGRALS OVER THE UPPER HALF SPACE (SPECULAR REFLECTION) 

Integrals over the upper half space (specular reflection) are the 
same as those over the lower half space except for the following. 
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